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ABSTRAGT 

We study the Kennicutt -Sclimidt st ar formation law and efficiency in the gaseous disk of the isolated 
galaxy GIG 96 (NGG 864: lTullyiri988| ). with special emphasis on its unusually large atomic gas (H I) 
disk (rHi/r25 = 3.5, r25 = 1.'85). We present deep GALEX near and far ultraviolet observations, 
used as a recent star formation tracer, and we compare them with new, high resolution (16", or 
1.6 kpc) VLA H I observations. The UV and H I maps show good spatial correlation outside the inner 
1', where the H I phase dominates over H2. Star-forming regions in the extended gaseous disk are 
mainly located along the enhanced H I emission within two (relatively) symmetric giant gaseous spiral 
arm-like features, which emulate a H I pseudo-ring at a r ~ 3'. Inside such structure, two smaller 
gaseous spiral arms extend from the NE and SW of the optical disk and connect to the previously 
mentioned H I pseudo-ring. Interestingly, we find that the (atomic) Kennicutt- Schmidt power law 
index systematically decreases with radius, from N ~ 3.0 ± 0.3 in the inner disk {Q'.S - 1.'7) to N = 1.6 
± 0.5 in the outskirts of the gaseous disk (3f3 - 4f2). Although the star formation efficiency (SEE), 
the star formation rate per unit of gas, decreases with radius where the H I component dominates 
as is common in galaxies, we find that there is a break of the correlation at r = 1.5 r25. At radii 
1.5 r25 < r < 3.5 7^25, mostly within the HI pseudo-ring structure, there exist regions whose SEE 
remains nearly constant, SEE ~ 10~^^ yr^^. We discuss about possible mechanisms that might be 
triggering the star formation in the outskirts of this galaxy, and wc suggest that the constant SEE 
for such large radii r > 2 r25 and at such low surface densities might be a common characteristic in 
extended UV disk galaxies. 

Subject headings: galaxies: spirals — galaxies: individual (NGG 864) — galaxies: structure — galaxies: 
ISM - stars: formation 
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1. INTRODUCTION 

Eor three decades it has been known that about 
one third of galaxies show unusu ally ex t ended 
atomic gas (H I) distribu tions (e.g., iBosmal 119811 : 
iHuchtmeier &: Richteil [l982| ) . It was not until recently, 
with the advent of the Ga laxy Evolution Explorer 
(GALEX) (iMartin et al.ll2005D . that we could observe 
easily star formatio n (SF) far b eyond the o ptical 
radius of galaxies (IThilker et all [20051 [20071 ^OOl 
IGil de Paz eralll200d 120071 : IBoissier ei~all[2007r i. The 
star-forming regions generally show a good spatial 
correlation with the H I component (the main phase of 
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neut ral gas at large ra dii) in such extended H I disks 
(e.g. iBigielet al.ll2010aD . 

These extended Ultraviolet (XUV) disks are character- 
ized by very blue UV-optical colors, they can reach up 
to four times the optical radius (as given by the /LtB=25 
mag/arcscc^ major axis isophote), and are mainly com- 
posed of low mass stellar clusters of 10'^ - 10^ M©. 
Recent studies have shown the presence of moderate 
amounts of dust and oxygen abundan c es of Zp^/lO 
(IGil de Paz et al.ll2007t iDong et al.ll2008t iBresolin et al.l 
|2009[) . iThilker et all (|2007D recently classified the XUV 
disks into Type 1 (p atchy UV emiss i on beyond the canon- 
ical SF threshold; iToomrd 119641: IMartin fc KennicuttI 
120011) and Type 2 (blue, large compared with the size 
of the galaxy in the near-infrared disks inside the SF 
threshold) and found XUV disks to be relatively frequent, 
20% and 10% of the overall disk-galaxy population, for 
respectively Types 1 and 2. 

Gomparing UV and H I in such extended disks allows 
the study of the Kennicutt-Schmidt (KS) law, or star 
formatio n rate (SFR) versus gas surface densities Tisfr 
vs Tsgas (|Schmidtlll9"5l lKennicuttlll998D . in the extreme 
low-density and often low-metallicity environments of the 
outskirts of galaxy disks, which can elucidate t he limit- 
ing c o nditions for gas cloud and star formation (jToomrd 
119641: IMartin fc KennicuttI [2001.1 . In addition, it is 
essential to study the relation between the availability 
of H I in the outer disk and the depletion time scales 
due to SF. Such an H I reservoir should play a role in 
the replenishment of the molecular gas content to enable 
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future SF. This process can thu s contribute to sustain 
SF over cosmological tim es (e.g. iKennicutt et all 119941 : 
iBauermeister et alll201ClD . 

Another key question to be addressed is what are 
the mechanisms for triggering SF in these fow-density 
gaseous environments. Several mechanisms have been 
proposed, both of external and internal origin, includ- 
ing gravit ational instabilities indu ced by galaxy interac- 
tions (e.g. iGil de Paz et al.|[2005[ ) or the i mpact of dark 
matter subhalos or extragalactic clouds (jTenorio-Tagld 
[1981: Bekki & Chiba 2005), turbul ence compression, su- 
perno vae, or gaseous s piral arms (lElrnegreen fc Hunten 
I2OO60 . In this regard. iThilker et al.l (|2007D found that 



75% of the Type-1 XUV-disk objects show some kind 
of optical morphology peculiarity or H I evidence for in- 
teraction/merger or a minor external per turbation, such 
as primordial clouds (jThilker et al.l 120091 ). For instance, 
NGC 4262, a member of the Virgo cluster, exhibits an 
extended H I and UV ring-like feature probably formed 
a s a result of a past m ajor interaction episode (?). 

IThilker et al.l ()2007() noted that isolated galaxies might 
also host a Type-1 XUV-disk. Although studying such 
isolated systems can provide provide information about 
other mechanisms producing extended disks, still no 
detailed study has been performed so far. An excel- 
lent t arget sample is the Catalog of Isolated Galaxies 
(GIG, lKarachentseval[T973[) . The isolation of the galax- 
ies in the catalog has been revised within the AMIGA 
project (Analysis of the inters tellar Medium of Isolated 
GAlaxies, http://arrLiga.iaa .es, IVerdes-Montenegro et al.l 
12001 iVerlev et al.l l2007allH ). The imposed isolation 
criterion assures that the galaxies have not suffered 
any m ajor interactions in a tirne sca le of the order of 
3 Gyr (jVerdes- Montenegro et al.l[2005[ ). We search ed the 
GAL EX Nearby Galaxy Survey (NGS, [Gil dc Pa z et al.l 
\m)% . Deep (DIS), Medium (MIS), and All Sky Surveys 
(AIS) for galaxies in the AMIGA sample. Out of the 24 
galaxies in the AMIGA sample having MIS, DIS and/or 
NGS data, none showed evidence of XUV emission. We 
found a total of 126 objects in AIS. Out of these, we 
identified seven XUV-disk galaxies, four of them of Type 
1 (GIG 96, GIG 103, GIG 712 and GIG 812). 

In this paper we examin e the properties of the ex- 
tended disk in the SAB(rs)c (jde Vaucouleurs et al.lll99i|) 
galaxy GIG 96 (NGC 864). GIG 96 has a strong and 
slightly curved bar, with a bar length est imated by 
Fourier analysis in H -band of 26" (jEskridge e t al. 2002; 
iSpeltincx et al.l[2"008[ ). It is locat ed nearby at a dista nce 
of 17.3 Mpc {V = 1561.6 km s-\ lEsnada et~an[2005l us- 
ing a IIo=75 km s~^ Mpc^^), which yields a linear scale 
1" ~ 100 pc. We adopted for GIG 96 the major and 
minor axis optical diameters D25 = 3.'7 and d25 = 2.'6 
at the isophotal level 25 mag/ arcsec^ in the B-ban d (ho- 
mogenized value from LEDA. [Paturel et al.ll2003l ). Note 
that the diameters in NED (NASA/IPAC Extragalactic 
Database) are larger, D25 = 4.'7 and d25 =3.^5, as in th e 
Uppsala General Catalogue of Galaxies (jNilsonI I1973D . 
However the homogenized diameters given by LEDA are 
co nsistent with r e cent v alues found in the literature, such 
as iPaturel et al.l (|2000( ) . We assume throughout the pa- 
per a position angle P.A. = 23°± 3° and an inclination i 
= 43°± 2° , values derived fro m modelling of the H I data 
cube in lEspada et al.l ()2005() . These are also consistent 



to the homogenized values presented in LEDA, P.A. = 
24.0° and i = 46.7°. 

VLA H I maps with resolutions of abou t ^ 45" reveale d 
an unusually extended H I disk (Esp ada et al.l l2005f) . 
Besides the optical revision of the isolation (jVerlev et al.l 
l2007a[ ). GIG 96 is also found to be isolated from any 
relatively H I-rich galaxy within the VLA primary beam 
(^^30'). Only a diminutive galaxy with an H I mass Mhi 
~ 5 X 10^ M0 located at 80 kpc from GIG 96 is its 
nearest companion. No tidal tails or other signs of inter- 
action were found to the reached sensitivity. Largest H I 
column densities were found within a ring feature with 
major and m i nor ax is of 1.'5 x IfO. More interestingly, 
lEspada et al.l (|2005[ ) found a pseudo-ring feature with an 
approximate size of 8.'2 x 4.'6. The galaxy H I disk has 
a kinematic asymmetry. The rotation curves in the ap- 
proaching and receding sides differ, mostly as a result 
of a kinematically detached clump which is evident as a 
secondary peak i n the position veloc ity cut at a radii r 
- 6.'5 to the SW (jEspada et al.l [20051 ). 

The main goal of this paper is to shed light into the 
local KS laws and Star Formation Efficiency (SEE), the 
Star Formation Rate (SFR) per unit of gas, of the ex- 
tended gaseous component in this isolated object, using 
new deep GALEX observations and high resolution {^ 
15") VLA H I observations (§[2]). In §[3|we describe our 
H I and UV maps, compare the spatial location of UV 
relative to H I, and study the KS law locally for different 
radius ranges, as well as the local SEE as a function of 
radius. In § |4| we compare the observed T,gas and rota- 
tion curve with that obtained from a model fitting using 
optical, UV and IR photometric radial profiles. Finally 
wc provide a discussion on the main findings and our 
conclusions in § [Sj 

2. OBSERVATIONS AND DATA REDUCTION 
2.1. VLA HI Data 

H I observations were carried out with the NRAO 
VLA[3 in its C configuration on July 23 2005, and 
were combined to th e D configuration data published by 
lEspada et al.l (|2005l) . We used the same correlator set- 
ting as for the D configuration data, which provided a 
velocity resolution of 48.8 kHz (10.4 km s^^) for the 64 
individual channels after Banning smoothing. The edit- 
ing and calibration of the data were done with AlPqj. 
The data were imaged using IMAGR by combining both 
datasets. For the cleaning process we applied a com- 
promise between natural and uniform weighting (weight- 
ing parameter ROBUST = 0). The average of the line- 
free channels has been subtracted from all the individual 
channels. 

The final rms noise level achieved after 4 (D- 
configuration) plus 7 hours (G-configuration) per chan- 
nel is (7 = 0.4mJy beam"^, with a synthesized beam of 
^ 16"9 X 15'.'6 (or 1.7 x 1.6 kpc), and a position an- 
gle P.A. = -30.10°. We obtained the integrated intensity 
(moment 0) map from the 27 channels with H I emission, 
to which we applied a primary beam correction (HPBW 
~ 30' for the VLA antennas) . The moment was masked 

^ The National Radio Astronomy Observatory is a facility of 
the National Science Foundation operated under cooperative agree- 
ment by Associated Universities, Inc. 

^ http://www.aips.nrao.edu/cook.html 
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using a 3cr (=1.2 mjy bcani^^) clipping for each single 
channel. The la noise level of the moment map 
is 0.02 Jy beam^^ km s^^. calculated as the addition in 
quadrature of the channel's la levels. 

The excellent agreement between the single-dish to- 
tal flux, S = 102.96 Jy k m s~^ using Green Bank 43m 
telescope (HPBW = 22') (jHavnes et all 1 19980 . and that 
obtained with the VLA C-l-D configuration data, S = 
100.3 Jy km s~^, implies that we recover most of the 
flux. Thus an extended gas component is not expected 
to be missing in our maps. 

2.2. GALEX NUV and FUV Data 

The GALEX observatory provides both a near UV 
(NUV, A = 1771 - 2831 A) and a far UV (FUV, A = 
1344 -1786 A) broadband filter, with angular resolutions 
(FWHM) of 5'.'3 and 4'.'2, respectively. The field of view 
of the instrument is about 1.25°, much larger than the 
VLA primary beam. CIG 96 was observed in NUV for 
6.385 ks and FUV for 1.648 ks in our GALEX Guest In- 
vestigator program 065 in Cycle 5. We also combined 
these data with previous observations from the archive. 

Foreground stars are prominent in the NUV images 
and we ide ntified and flagged the m easily via their UV 
colors (see iGil de Paz et al.l l2007f ) . We computed and 
subtracted the background from both the FUV and NUV 
maps. We corr ected fluxes by adopt ing Galactic extinc- 
tion as given bv lSchlegel et al.l ()1998n . namely E(B-V) = 
0.048. 

Besides correcting for the foreground galactic extinc- 
tion, we must also account for radial variations in the 
internal extinction. The total IR (TIR) to UV ratio, 
TIR/UV, is a robust tracer of internal extinction, de- 
pending weakly on th e dust-to-stars geometry or the ex- 
tinct ion law (see, e.g. iWitt fc GordonI 120001 : iBuat et al.l 
120051 and references therein). In order to compute 
the TIR radial profiles, Spitzer images at 24, 70 and 
160 ^m would be required, but they are not available 
for this galaxy. Therefore, we resorted to the FUV- 
NUV colour (i.e., the UV slope) as an indirect tracer 
of dust extinction. Instead of using th e classical cal- 
ibrations derived fo r starburst galaxies (iCalzett i et al.l 
[1994; He ckman et aL 1995: Meurer et all [l 9991. here w e 
made use of the recipes of iMuhoz-Mateos et al.l (|2009D , 
which were calibrated on a sample of normal, nearby 
galaxies. From the FUV-NUV colour profile we in- 
fer the TIR/FUV and TIR/NUV ratios at each radius. 
These ratios are then converted into extinction radial 
profiles in the FUV and NUV bands indep endently. This 
was d one by means of the prescriptions of iCortese et al.l 
(|2008f ). which take into account the varying contribution 
of young and old stars to the dust heating. The error in 
the extinction is dominated by the scatter of ratios when 
we cahbrate TIR/FUV and TIR/NUV as a function of 
colour FUV-NUV, and is estimated to be of the order of 
0.4 mag. 

3. RESULTS 

3.1. HI Component 

We present in Fig. [T] (lower left panel, green) the H I 
integrated density distribution. The H I extent is Djji 
= 12f9 at a level 0.7 Mq pc'^ {la noise level, § [2T|) . 
which results in an unusually large H I to optical extent 



ratio of DH1/D25 ~ 3.5. A similar ratio is that of M83, 
whose X UV disk properties have been previously studi ed 
in detail (|Thilker et al.l[2005L [2007t IBigiel et al.(l20T0bh . 

The VLA C-l-D configuration H I map reveals a large 
variety of structures that could not be clearly distin- 
guished in the VLA D-configuration map due to the lower 
angular resolution. The inner ring, with projected diam- 
eters of 1.'5 X l.'O, exhibits enhanced H I emission along 
the stellar spiral arms. The outer pseudo-ring, with an 
extent of 8.'2 x 4f6, seems to be connected to the in- 
ner ring to the NE and SW of the inner disk. Overall 
the shape of the outer pseudo-ring resembles two giant 
spiral arms surrounding the optical disk. 

Although the shape of the brighter features is quite 
symmetric, there is also a faint diffuse H I emission that 
is asymmetric. The kincmatical ly detached c o mpon ent 
in the SW that was reported in lEspada et al.l (|2005| ) at 
about a radius 6' is seen here as a clumpy and extended 
component. A more detailed analysis of the whole H I 
data-cube will be presented in a forthcoming paper. 

3.2. FUV and NUV Emission 

We present in Fig. [T] the NUV map (upper left panel, 
red) and FUV map (upper right panel, blue) with the 
same field of view than that of the H I map. Both NUV 
and FUV maps show an almost identical distribution. 
However, while the radially averaged FUV-NUV color 
^ 0.5 in the internal parts of the disk (in units of AB 
magnitudes per square arcsec), it reaches FUV-NUV = 
0.1 - 0.2 in the most external regions. 

First, there are two nearly symmetric arms, one start- 
ing at the E of the bar and extending towards the N, 
and another one starting at the W of the bar and ex- 
tending toward the south. The SW arm's appearance in 
both NUV and FUV emission i s clumpy. This i s sim ilar 
appearance to that o bserved bvlEskridge et all ()2002f ) in 
the optical and NIR. lEskridge et al.l claim that the spi- 
ral pattern is flocculent in the outer parts of the inner 
optical/NIR disk (r < 1.'5). These fiocculent regions are 
related to the connections between the outer H I pseudo- 
ring and the inner H I ring to the NE and SW at the P.A. 
of the inner disk. A central concentration is found unre- 
solved, which is like ly connected to a nuclear starburst 
(jMartini et al.|[2003[ ). 

Second, a more diffuse and extended component of UV 
emission is seen in the UV maps. Two symmetric spiral 
arms are apparent, extending from the NE and SW of the 
optical disk. These spiral arms join an outer pseudo-ring 
structure to the E and W starting at about r = 2'. The 
pseudo-ring structure is the same as that seen in the H I 
outer disk, although the star-forming regions seem to be 
more confined to regions with high H I surface density. 
Other filamentary and more chaotic structures are seen 
out to almost r = & , especially to the SW and to a minor 
extent in the NE. 

A composite image of the H I and NUV, FUV maps 
convolved to the resolution of the former is presented in 
Fig.[l](lower right panel). An scheme (over the composite 
image) of main features in the outskirts of CIG 96 is 
shown in the upper panel of Fig. [2j spiral arms (blue 
cross signs) extending from the N to the W and from the 
S to the E, connecting the pseudo-ring structure (red 
cross signs). The circle signs indicate the SE spiral arm 
and pseudo-ring (SE side too) rotated by 180°, which 
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shows that these features are nearly symmetrical. 

The lower panel of Fig. [2] shows that the NUV (colour 
pixel map) and H I maps (contour map) have a high de- 
gree of correlation. This tight spatial correlation of H I 
and UV emission remains far into the outer pseudo-ring, 
and a diffuse component all over the H I envelope, as it 
can be seen in the last plotted contour. 

3.3. Star Formation Law 

We obtained the KS law (jSchmidtl 119591: iKennicuttl 
|1998[) pixel-by-pixel by comparing the gas sur face den- 
sity (Yigas) map obtained from our H I data (§ 12. ip and 
the SFR surface density (SgpR) from the extinction- 
corrected (radially) NUV data (§ 12. 2p . We converted 
the pixel size of both maps to 16" (as the H I beam size). 

In the outskirts of disks the gaseous component is usu- 
ally dominated by H I rather than by molecular gas, 
so the H I content is a good estimate of the total gas 
content, Sgas — Sjj j, especially below the saturation 
limit between atomic and molecular gas at 'S,gas — 
12 M ^ pc-^ (iMartin fc KennicuttI [200lt IWong fc Blitz! 
[200llBigiel et al.ll2008l: IVerlev et al.ll2009L 120101) . 

In order to calculate the Egas we took into account 
the inclination of the galaxy and we multiplied by 1.36 
to include the He lium contribution ([Walter et al.l [20081 : 
iLerovet al.]l200l : 

Sg„^[M0pc-2] = 33.6 X EiSi[Jybcam-i]AV[km s'^] 

where Si is the flux density for each channel and AV is 
the channel velocity width. The la noise level for Egas 
= 0.7 Mopc~^, using the sensitivity limit reached in 
our H I moment map (§ HTJ. 

As for the 'Ssfr map we used the corrected NUV data 
(§ 12. 2p . In the literature, FUV is more widely used in 
order to derive SFR as it involves emission from young 
stars born in the last 100 Myr, more similar to the Ha 
< 10 Myr) than to the NUV emission (< 300 Myr) 
Verlev et al.ll2009l ). However, by using NUV, we achieve 
a gain in sensitivity and NUV emission is not essentially 
different to FUV to our resolution. This approach can 
also be useful for other sources with more sensitive NUV 
data than FUV. As for the 'hidden star formation', i.e. 
the SF that we cannot correct from the UV color be- 
cause it is completely obscured by dust in both bands, 
is ex pected to be negligible in th e outer parts of disks, 
since iPrescott fc Kennicuta ([2008f ) showed that there are 
no such regions at rad ii r > 0.6 ^25 for SINGS galaxies 
([Kennicutt et al.ll2003[) . 

We convolved the NUV map to ma t ch th e resolution 
of the H I map. Following [KennicuttI ([1998f ). the Esfr 
was calculated as: 



^SFB.[Mq yr-l kpC^^] = IQ-O-^XMNUvlABmag arcsoc ]+7.41 

assuming solar metallicity, a lSalpeten ([1955| ) initial mass 
function (IMF) and that the SFR has remained constant 
over the last few 10^ yr, which is the typical lifetime 
of stars dominating the UV emission. In this equation 
IJ'NUV is the surface brightness (in AB magnitude units) 
of NUV emission. To convert to the SFR derived using 
iKroupal (l2001[ )'s IMF one should multiply our SFR by a 
factor 0.83. The Ict sensitivity limit is ^sfr = 9.2 x 



10~^ M0 yr"-'^ kpc~^. In order to avoid projection effects, 
we have also deprojectcd both the "Ssfr and the Eg^s 
maps. 

Fig. [3] shows the correlation between the ^sfr de- 
rived from both the NUV and FUV maps pixels (16"size) 
within the inner 400". The J^sfr obtained from the FUV 
map has been calculated in a similar manner as the T,sfr 
from the NUV map. The error bars presented in Fig. [3] 
represent the formal uncertainties obtained from the sen- 
sitivity limit as well as from the uncertainties associated 
to the extinction correction. A bisector fit to all the data 
points yields a slope equal to unity and an intercept close 
to (1.02 ± 0.03 and 0.15 ± 0.01, respectively). We 
estimated the uncertainties using different random rep- 
resentations of the data taking into account the uncer- 
tainties (bootstrapping). We obtained similar uncertain- 
ties for the formal slope and intercept, 0.01 and 0.02, 
respectively. Finally, note that calculations of "Ssfr us- 
ing al ternate SFR tracers give differences of the order of 
50% ([Bigiel et al.ll2010al: iLerov et al.l [20081: IVerlev et al.l 
|2009() . We do not consider this uncertainty in our error 
estimation. 

In Fig. m we show the T,sfr vs Eg^s plot for GIG 96. 
From the figure we can discern the saturation limit be- 
tween atomic and molecular gas found for other galaxies 
at Ega, ~ 12 Mq pc-2. 

In order to inspect whether the local KS star for- 
mation law changes as a function of radius, we cal- 
culated a linear fit for different annuli centered at 
a(2000) = 02^^ 15" 27"64 (5(20 00) = 06° 00' 09" 1 
([Leon fc Verdes-Montenegroll2003D . Each colour in the 
data points of the ^sfr. versus Eg^s plot in the lower 
panel of Fig. [4| represents a different annulus zone over 
the deprojectcd H I map, as depicted in the upper panel 

of Fig. m ' 

First, in the innermost region (r < 0f8, or 4.8 kpc; 
in dark blue in the left upper panel of Fig(3|), the cor- 
relation between T,sfr. versus T,gas is biased since at 
such radii the contribution of gas in molecular phase 
is presumably larger than that in atomic phase. For 
instance, Mh2 = 2 x 10* Mg using CO(l-O) observa- 
tion at the Swedish-ESO 15m Submillimeter Telescope 
(SEST), which is characterized by a HPBW = 42" (ra - 
dius - 2 kpc) ([Elfhag et al.|[l99l[Lismfeld et al.|[201l[) . 
If the molecular gas is distributed within a uniform dis- 
tribution then the corresponding molecular gas surface 
density would be E//2 ^ 11-5 M© pc~^. The Y,h2 is 
likely larger since enhanced emission is expected toward 
the nucleus, inner spiral arms and along the bar. Should 
we include the molecular gas contribution, the correla- 
tion T,SFR ~ Egos would likely continue above Eg^s > 12 
Mq pc~^ in a standard manner, with a slo pe N ~ 1 - 
1.4 a s it is observed in other galaxies (e.g. iBigiel et al.l 
[2OO8I) . 

We find that the power index (bisector fit) of the KS 
law systematically changes with radius from N = 3.0 ± 
0.3 in the inner disk (O'.S - 1.'7, or 4.8 ~ 10.2 kpc) to 
N — 1.6 ± 0.5 in the pseudo-ring feature (3.'3 - 4f2, or 
19.8 - 25.2 kpc). Further out the correlations are bi- 
ased since a considerably large amount of data points 
fall below the sensitivity limit. We detailed in Table [Tj 
the parameters for the fit with T,gas as independent vari- 
able (OLS(Y — X)) as well as the least square bisector 
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fit, and for annuH from 0.'8 - 4f2 in bins of 0f8 (4.8 - 
25.2 kpc in bins of 4.8 kpc). The estimated errors of the 
fit parameters have been calculated via bootstrapping, in 
the same manner as in the T,sfr comparison from NUV 
and FUV. Note that in general these param eters are not 
strong ly affected due to spatial resolution. iBigiel et al.l 
(|2008[ ) shows that the power law parameters vary only 
weakly with changing spatial resolution. 

Besides the uncertainty that results from the assump- 
tion "Egas — Sjj J, note that there is another source 
of uncertainty regarding a possible metallicity gradient 
that might affect the calculated Ssfj?, (in the sense that 
metallicity is lower in the outskirts). Assuming a con- 
stant calibration for a given metallicity, it can be inferred 
that regions with a lower metallicity wo uld have a lower 
S.SFR for a fixed UV surface brightness (jLeitherer et al.l 

Mm . 

3.4. Star Formation Efficiency 

The SFE, the SF per unit gas mass, decreases with 
radius from values of 10"^*^ yr""'^ in the inner ring to ~ 
10^11 yr~i in the outskirts. In Fig. [5] we show the SFE 
as a function of radius (normalized to r25—D25/2). An 
increase in the observed scatter is present as a function 
of radius, from ^ 0.2 dex to ~ 1 dex. 

For comparison we have plotted in Fig. [Slthe derived 
fit for the SFE as a function of radius frorri iLerov et al.l 
11)08): 

SFE = 4.3 X 10"^" yr~\ ff r < 0.4 ras 

SFE = 2.2 X 10"^ exp(-r/(0.25r25)) yr"\ff r > 0.4 T25 

Note that their calibration is valid up to 1.2 r25 
(|Lerov et al.ll2008| ). 

The agreement between this fit from iLerov et al.l 
(|2008[ ) and our SFE radial profile is reasonably good at 
0.4 r25 < r- < 1.5 r25. Although our fit is slightly offset 
towards higher values, the slope is essentially the same 
in this radius ran ge. The fit to our data and that of 
ILerov etlUI (|2008| ) are shown in Fig. [5] as d ashed lines. 

From r > 1.5 r25 the fit bv ILerov et "all (|2008| ) would 
start to deviate considerably from our data, in the 
same manner as po inted out in other spiral galaxies by 
IBigiel et al.l ()2010a[) . It seems to be more constant, at 
least for those regions where Ssfr and E^as sue above 
the noise level (blue data points). For such regions the 
SFE ~ 10"". The fit to our data is shown in Fig. [5] 
as a green dashed line. It fits relatively well (rms ^ 0.5 
dex) over a radius range 1.5 r25 < r < 3.5 r25. Note that 
the scatter is large enough to consider the fit parameters 
just as estimates. This region corresponds to the pseudo- 
ring feature. Many of the detected regions in the radius 
range 1.7 r25 < r < 3.5 r25 are with in the kinematically 
detached clump located to the SW (jEspada et al.1120051 ) 
(green circles). 

Red triangle signs in Fig. [5] correspond to upper limits 
in TisFR, although detected in Sgas, and thus upper lim- 
its in SFE. These data points would probably increase 
the scatter of SFE at large radii. 

4. MODEL OF THE MULTI- WAVELENGTH LIGHT 
DISTRIBUTION OF CIG 96 

We inspect in this section whether the photometric 
properties, including the XUV emission, is common of 



a spiral galaxy disk (with a given spin and mass), or 
whether on the contrary it is a peculiar object not only 
for its photometric properties but also for its SF history. 

4.1. Model 

We derived photometric radial profiles for CIG 96 not 
only in the GALEX FUV and NUV bands, but also in 
the Sloan D igital Sky Survey ( SPSS) ugriz bands and 
the Sp itzer (jWerner et al.l |2004) IRAC 3.6 and 4.5 urn 
bands (jFazio et al.l(200l 0. 

These multi- wavelen gth radial pr ofiles were fitted with 
the model s of Hp issier fc PrantzosI (|l999. 2000) with the 
IMF from iKroupal (|2001[ ). which are able to predict the 
radial variation of several galactic properties as a funcion 
of two parameters: the spin A (non-dimensional specific 
angular momentum) and the asymptotic velocity of the 
rotation curve Vc- Galaxies are modeled as a set of in- 
dependently evolving rings, with the gas infall time-scale 
depending on both the total galaxy mass and the local 
mass surface density. Gas is then converted into stars fol- 
lowing a Kennicutt-Schmidt law, modulated by a dynam- 
ical ter m that mimics the per iodic passage of spiral arms 
(see also lBoissier et al.ll2003a[) . Newly-bor n sta rs follow a 
user-specified IMF (here we rely on that of lKrou pa 2001). 
Stars of different masses die at different rates, enriching 
the ISM with metals that are incorporated in subsequent 
generations of stars. The lifetimes, yields, evolutionary 
paths and spectra of stars at different radii are computed 
as a function of the local metallicity. An initial model 
was first calibrated to reproduce observables in the Milky 
Way (Boissier & Prantzos 1 999i). This model was the n 
generalized to other galaxies (jBoissier fc Prantz os 2000) , 
using the sca ling laws derived from the A-Cold Dark Mat- 
ter scenario (jMo et al.1119981 ). 

The GALEX, SDSS and Spitzer photometric profiles 
were simultaneously fitted with these disk evolution mod- 
els. The internal extinction profiles in the optical and 
near-IR bands were derived from the UV ones using a 
MW extinction law and a sandwich dust-to-stars geome- 
try (further details on the fitting procedure can be found 
in lMunoz-Mateos et al.|[2"011l ). By minimizing the x^ be- 
tween the observed and predicted photometric profiles, 
the best-fitting values were found to be A = 0.048(|q °2^'^ 
and Vc = 16llg^km s^^. Using other models character- 



ized by an IMF with less massive stars (i.e. iKroupa et al.l 
Il993f) or excluding the XUV disk further out than 150" 
does not successfully reproduce both photometric profiles 
and the observed H I ro tation curve as the ones used here 
do (jEspadaet al.ll2005D . 

We plot in each panel of Fig.|6]the de-projected surface 
density profile for each band considered, as well as the 
fits from the model. The agreement is good for all com- 
ponents. For UV bands, the azimuthally averaged AB 
magnitudes seem to be slightly underestimated in the r 
= 10"- 15"region, although in general they agree well. 

The KS law was expressed in §[3?3]in the form Y^sfr oc 
E^jj, with N varying for different annuli. On the other 
hand, the SFR in the model is expressed as Tisfr. oc 
^g'as V[R)/R. To check the consistency between both 

•^ We used the Basic Calibrated Data publicy available in the 
Spitzer archive. These images were obtained as part of the Spitzer 
Survey of Stellar Structure in Galaxies (P.I. K. Sheth). 
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definitions, in Fig. 2] (bottom) we plot the temporal evo- 
lution from (0 to 13.5 Gyr) of Esfr and Sgas using the 
model, and for different radii, 5.5 (blue), 10.9 (green) 
and 20 kpc (red). No obvious offset is found, which indi- 
cates that this galaxy is in general following the SF law 
implemented in the models, which relies on average main 
properties of galaxies. 

Note that the models used here are not necessarily 
good at reproducing the properties of galaxy disks out 
to these large radii. The model itself does not include an 
imposed threshold, and although it is an extrapolation, it 
should work reasonably well for low density regions and 
outskirts as long as there is no strong threshold effects 
and the mode of star formation does not change. In fact 
these models have been tested ag ainst observations of 
Low Surface Brightness galaxies in lBoissier et al.ll2003bl 

4.2. Predicted Rotation Curve and Gas Surface Density 

The model provides an estimation of the rotation curve 
and 'Sgas as a function of radius. We compare the ob- 
served gas surface density with that predicted by the 
model in the upper panel of Fig. [T) A small deviation is 
found at low radii, where most of the gaseous component 
is in molecular phase, and at large radii from r > 15 kpc 
(> 3') where the model underestimates the observed sur- 
face density as a result of the external pseudo-ring and 
the extended UV emission in general. 

The latter can be probably explained as a result of 
a past recent event that re-distributed part of the gas 
on time scale shorter than that probed by the differ- 
ent bands, because they average over the lifetime of the 
galaxy. UV emission, which as H I should be sensitive 
to any perturbation, agrees relatively well to the model. 
This suggests that a recent event could have occurred, 
such as accretion, in order for the H I to show this devi- 
ation with respect to the model. 

In the lower panel of Fig. [7] it is also seen that 
the agreement betw een the observed H I rotation curve 
(jEspada et al.l [20051 ) and that predicted by the model is 
excellent. Note that the H I observed rotation curve is 
not used during the fit. 

5. DISCUSSION AND CONCLUSIONS 

The outer regions of disk galaxies are essential as test- 
ing sites for models of disk assembly and evolution. 
XUV-disk galaxies are ideal laboratories to study low 
SFR and low gas surface densities conditions. SF in 
the extreme (low-density, low-metallicity) environment 
of the outskirts of galaxy disks can elucidate the limiting 
conditions for gas cloud formation and subsequent SF. 

GIG 96 is a unique object because it resides in isola- 
tion (> 3 Gyr since last major interaction) and has an 
extended star forming disk, even though XUV disks are 
usually ascribed to major interaction events. The ob- 
served gas surface density from 2'.5 to 4.'7 (15 kpc to 
28 kpc) is larger than that expected from our model, 
even though the kinematics is in average well repro- 
duced. This suggests that this extra gas and SF could 
be an anomaly in a (until-recently) well-behaved disk 
with a (non-dimensional) specific angular momentum 
A=0.048^0Qi[2 and an asymptotic velocity of the rota- 
tion curve Vc = 161]']g^ km s~^. 

Our sensitive and high resolution data allowed us to 



derive the power law index of the KS law locally for dif- 
ferent annuli, even at large radii. We find that the slope 
decreases with radius, being N ~ 3.0 ± 0.3 in the inner 
disk (O.'S - lf7) and decreasing to N = 1.60 ± 0.5 in the 
outskirts of the gaseous disk (3f3 - 4'2). The KS law 
index tends to unity as radius increases. This might re- 
veal a privileged relation between H I and SF where the 
amount of molecular gas content traced by CO is small. 
H I could thus be a good tracer of star forming gas in 
this low density low metallicity conditions. More sensi- 
tive observations of both H I and UV are necessary to 
reveal what is the limit for the KS index in extended 
disks. 

We found that up to the pseudo-ring feature the SFE 
radial profile is typical of other disk galaxies, with SFE 
decreasing with radius. Overall, the SFE spans typical 
values from '^ 10~^ yr~^ to 10~^^ yr~^. Interestingly, 
the outer part of the H I disk, the pseudo-ring feature, 
shows a break in the decreasing t r end o f the SFE pre- 
sented for example by iLerov etall (|2008f ). The SFE for 
regions outside the inner disk (r > 1.5 r25) is SFE ~ 
10~^^ yr~^ (§ 13. 3p . There are regions that share this 
consta nt value up to 3.5r2fl , especially the southern part 
which lEspada et all (|2005i ) showed that it is actually a 
kinematically detached region to the rest of the disk. Re- 
gions with no relevant T^sfr but above the noise level 
of T,gas suggest very low SFE, of the order of about 
10"^^ yi'~^- The depletion time scales involved are quite 
large and thus the gas is very likely a reservoir for future 
SF. Even for the pseudo-ring feature that has a larger 
SFE for its radius, if 1% of the gas turns into stars in 
10^ yrs, the time scale for consuming all the gas would 
be large, of the order of ~ 10 Gyr. 

Are these deviations in SFE at large radii the rule more 
than the exception in other XUV disks? There are two 
pieces of evidence that seems to show that these devia- 
tions from normal values m ay be commonly pr esent in 
other galaxies. In Fig. 1 of ILerov et al.l (|2008f l. where 
SFE is plotted vs galactocentric radius, the average of 
SFE for all the data points is larger than their fit for r/r25 
> 1.0. Unfortunately their calibration limit is r/r25 = 
1.2. Probably the most relevant case to compare the SFE 
at large radii is that of M83. M83's extended H I disk 
reaches, as in GIG 96, a radius r ^ 3.5 r^F, , and it has 
been w idely studied in the UV and H I by iBigiel et"al] 
(|2010bD . The SFE in both cases level off at r ~ 1.5 r 25 
with a value of 10~^^ yr~^. Note that for M83 this cor- 
responds to th e depletion time Td^ j, = 10^ Gyr in Fig. 4 
right panel of IBigiel et al.l l2010bf ) , also with a scatter 
0.5 - 1 dex. These evidences suggest that we might have 
higher SFE in extended disks with respect to that extrap- 
olated from inner radii. In order to check the generality 
of this result, deeper H I and UV observations of more 
XUV disks are needed. 

Next question is what mechanism is enhancing the 
H I densities and triggering the SF in the outskirts of 
these galaxies? The isolation of GIG 96 since about 
~ 3 Gyr ago allow us to discard the major interaction 
event scenario. Thus, if g ravitational instabilitie s in- 
duced by interactions (e.g. iGil de Paz et al.l 120051 ) are 
producing the observed spiral arms and pseudo-ring fea- 
ture, then we can reject the major interaction event. It 
must be due to minor companions (hypothesis explored 
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bv lEspada at al.ll2005D. extragalactic clouds, or dark mat - 
ter sub-halo s (iTenorio-TaglelllQSlUBekki fc Chiba|[2006l ). 
lEspada at al.1 ()2005D suggastad the possibility of an ac- 
cretion of a Hl-rich companion to be the responsible of 
the SW kinematically detached clump and the overall 
outer ring-like structure. 

However, given the remarkable spiral symmetry found 
in the UV and H I maps (i.e. two symmetric spiral 
arms from the N and S of the optical disk, and con- 
necting the outer pseudo-ring feature to the E and W, 
respectively), one may think that it is an intrinsic prop- 
erty of the galaxy rather than induced by the environ- 
ment. The outer spiral arms could be due to an insta- 
bility or could be driven, either by a non -axisymmetric 
halo fe.g. iMachado fc Athanassoulall2010l ). or by a non- 
axisymmetric disk. We favor this last option, although 
it is difficult to prove, or disprove since the major axis 
of this oval would be more or less aligned with the kine- 
matical major axis. But the global morphology makes 
this option quite plausible. In particular, the arms em- 
anate symmetrically in the NW and SE, and then fall 
back at the other si de. This is a pa ttern expected in 
the manifold theory (iRomero-Gomez et al. 2006, 2003; 
lAthanassoula et al.ll2009aU E l2010[ )lfthe disk is oval and 
the manifolds allow the orbits to make a 270° turn, while 
meeting the opposite arm. 
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Fig. 1. — GALEX near UV (upper left, red) and far UV (upper right, blue) maps, as well as the VLA (C+D configuration) H I moment 
map of GIG 96 (lower left, green), and a composite of the three maps (lower right) where the resolution of the UV maps have been 
degraded to that of the H I. All the panels cover a field of view of 13' X 13'. The resolution of each map is shown in the lower right. 



TABLE 1 
Local Kennicutt-Schmidt law for different annuli 



OLS(Y— X) 



Bisector 



Zone Radius range (') / (kpc) Slope (N) Intercept Slope (N) Intercept 

i 0.8 - 1.7 (4.8 - 10.2) 1.7 ± 0.3 

2 1.7-2.5(10.2-15.0) 1.9 ± 0.3 

3 2.5-3.3(15.0-19.8) 1.4 ± 0.2 

4 3.3-4.2(19.8-25.2) 1.2 ifc 0.2 



-4.0 ± 0.3 3.0 ± 0.3 -5.4 ± 0.3 

-4.9± 0.2 2.18 ± 0.10 -5.07 ± 0.11 

-4.85± 0.14 1.7 ± 0.2 -5.1 ± 0.3 

-5.1± 0.10 1.6 ± 0.5 -5.3 ± 0.6 
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Fig. 2. — Left panel) Scheme (over the composite image) of the gaseous spiral arms (blue cross signs) extending from the N to the W and 
from the S to the E, and the pseudo-ring structure (red cross signs). The circle signs indicate the SE spiral arm and the SE pseudo-ring 
side, rotated by 180°. Right panel) H I map in contours over the NUV color (logarithmic) scale map. The H I contours are: 0.02, 0.08, 
0.16, 0.24 and 0.32 Jy beam"! km s"!, or Sjj j = 0.5, 1.9, 3.9, 5.8, 7.8 Mq pc'^. 
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Fig. 3. — Pixel-to-pixel (16") comparison of SFR obtained from both extinction corrected NUV and FUV maps. The (red) line indicates 
the bisector fit using all the data points above the sensitivity limits, with a slope equal to unity (1.02 ± 0.03) and an intercept close to 
(0.15 ± 0.01). The two dashed lines represent the nois e level limit for the SFR derived from NUV (horizontal) and FUV (vertical). Error 
bars have been calculated as explained in Section 13.31 
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4. — Upper panels) Deprojected H I map (left), with the colour scale showing different radial regions (right). Lower panel) SgpR 
gas (KS law) pixel to pixel for the different rings. Linear fits (bisector) to KS law for the different regions are also shown as dotted 
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Fig 
vs S 

lines (See Table[TJ. The diagonal (red and blue) dashed lines indicate constant SFE levels at 10"^'^ (bottom line), 10~^^, 10"^" 
and 10~* (upper line) yr~^. The vertical thick dashed line shows the threshold between molecular and atomic gas at Sgas = 12 Mq pc 
The vertical and horizontal thick dotted lines indicate the noise levels of Sgas (vertical) and EsFR (horizontal), at O.TMq pc~-^ and 1.9 X 
10~^Mq yr"^ kpc~^, respectively. The cross signs represent the data points in the kinematically detached region in the SW l lEspada et al.l 
[2005f) . 




Fig. 5. — Star Formation Efficiency (SFE) (in logarithmic scale) as a function of galactoccntric radius normalized by r25, up to r/r25 = 
3.5. Blue data points indicate regions whose Sgas and SgpR both are above their noise levels, as in Fig.|4] There is an apparent break in 
the SFE fit at around r/r25 = 1.5. The green dashed lines represents a fit to the detected data points: intercept = -8.2 ±0.02, slope = 
-1.69 ill 0.02 if r/r25 < 1.5 and intercept = -10.5 ±0.1, slope = -0.14 ± 0.05 if r/r25 > 1.5. The grey dashed lines show the fit obtained by 
ILerov et al.l (2008) for 18 nearby galaxies (see ^ 13.41 1. Below r/r25 <0.4 the gas surface density is likely dominated by molecular gas and 
the correlation to the fit does not hold. Red data points correspond to regions whose Sgas is above the noise level but Ssp R is an upper 
limit. Green circles represent the detected data poi nts in the kinematic ally detached region in the SW l|Espada et al.|[2005l) . The black 
triangle signs represent the median data points from lBigiel et aLl Il2010al ) for a sample of spiral galaxies. 
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Fig. 6. — Surface density radial profiles for the different bands: GALEX NUV and FUV, SDSS ugriz bands, and Spitzer 3.6 and 4.5 ^m), 
in bins of 6". The grey data points are corrected for Galactic extinction and the black data points for internal extinction as well. The 
latter are used for the fits using a model based on the IMF by IKroupal l|20011 ) (A = 0.048_qqj2 and Vc = 161_9 , as explained in §|4ll. 
The best radial profile fits are shown as (red) solid lines. The red band indicates an estimation of the error, and comprise all the fits within 
2 times the minimum x^ ■ The red dashed line indicates the radius we use to separate the region of the profile dominated by the bulge 
(left) and that by the disk (right). Only the latter region is used to perform the fit. A scale conversion for each band to Sspr in units of 
Mq yr~^ kpc~^ is also shown in the upper left panel. 
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Fig. 7. — Upper panel) Observed gas surface density (Egas) radial profile and fit using our model (§ ll} based on the IMF by lKroupal 
11200 II) . Lower panel) Comparison of tfie observed rotation curve (average southern and northern sides, [Espada et al.l 120051 ) and that 
predicted by our model. 



